Coherent excitation transferring via dark state in light-harvesting process 
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We study the light absorption and energy transferring in a donor-acceptor system with a bionic 
structure. In the optimal case with uniform couplings, it is found that the quantum dynamics of this 
seemingly complicated system is reduced as a three-level system of A-type. With this observation, 
we show that the dark state based electromagnetically-induced transparency (EIT) effect could 
enhance the energy transfer efficiency, through a quantum interference effect suppressing the excited 
population of the donors. We estimate the optimal parameters of the system to achieve the maximum 
output power. The splitting behavior of maximum power may be used to explain the phenomenon 
that the photosynthesis systems mainly absorb two colors of light. 

PACS numbers: 71.35.-y, 73.22.-f, 87.15.Mi 



To tackle the global problem of energy source [l[ peo- 
ple may learn lots from the natural process of light- 
harvesting in plants, algae and bacteria. The energy 
transfer mechanism of high efficiency in light-harvesting 
process would help to design the new generation of clean 
solar energy sources. Recently, the long-time coherent 
properties of excitation in light-harvesting systems have 
been observed in experiments. This coherence can be pre- 
served in these structures for a long time, even in room 
temperature. It seems that the high efficiency of light 
conversion is related to these coherent properties even 
with quantum natures. Therefore, the physical mecha- 
nism of photosynthesis assisted by (somehow quantum) 
coherence attracts much more attention from both ex- 
perimental [2] and theoretical aspects 

An inherent mechanism for high efficiency energy 
transferring may be due to the optimizing of spatial 
structure of light-harvesting systems. The X-ray analysis 
(lij has revealed some common elements shared by differ- 
ent light-harvesting complexes in nature, one of which is 
a ring structure with a centered reaction center [l7l [l8| . 
It is well-known that the nature selection rules always 
keep the most adaptable feature in the biological sys- 
tem for the present environment. Therefore, it is be- 
lieved that this kind of structure takes advantage in the 
light-harvesting process. The mechanism behind this op- 
timal structure may account for the high efficiency of 
the natural light-harvesting process. Thus the investi- 
gation of mechanism of similar system would be heuris- 
tic to design artificial light-harvesting systems with self- 
assembling molecular array or the quantum dot array in 
top-down-semiconductor fabrication in the future. 

In this letter, we will study a generic model, which is 
similar to the light harvesting complex of type I (LHC I) , 
a centralized acceptor surrounded by the coupled donors 
arranged in a ring. In natural light-harvesting process, 
the light-capture process and excitation transferring hap- 
pen simultaneously. To mimic the natural process, we 
include light capture process in the model by coupling 
the donors with photons in single mode. By showing the 
present model with homogeneous coupling could be re- 
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FIG. 1: (Color online) (a) The photon assisted donor- 
acceptor system with dissipations in both donors (rate k) and 
acceptor (rate T); (b) The reduced energy spectrum in single 
excitation subspace for its structure optimization 



duced into the well-studied three-level A system, we find 
that the overall transfer efficiency is insensitive to the 
decay of the donors when the eigen-frequency of the ac- 
ceptor is resonant with the light. This discovery implies 
that the dark state effect suppresses the excitation popu- 
lation on the noisy donors so that the energy transferring 
efficiency is dramatically improved. Otherwise, the loss 
of donor excitations will largely decrease the transferring 
efficiency to the acceptor. 

As illustrated in Fig Ufa), our model concerning light 
capture and coherent excitation transfer, is similar to the 
structure of LHC I. The one-dimensional circle array con- 
sists N donors which is also analogous to the ring struc- 
tures in the LHC II, such as B800 or B850 ring. Similar 
to the pigment molecules in natural LHC I ring, each 
donor can be modeled as a two-level system (TLS) |e,) 
and \gi) with energy level spacing (i = 1, N). TLS 
is a proper approximation for single excitation case in 
practice. In most natural_conditions, the LHCs always 
have only one excitation 



thus our discussion only 



concerns the zero and one-excitation subspaces. The ac- 
ceptor is placed at the center of ring. Since the hoping 
between the non- adjacent sites would be weak, we con- 
sider only the adjacent hopping of excitation on the ring 
with strength g. As we know, the visible light takes main 
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part of the energy of the solar spectrum and the corre- 
sponding wavelength is about 5 x 10 3 A. In the cell of 
photosynthesis bacteria, the B800 and B850 ring usually 
have the radius about 4 — 6nm 17]. Thus, the BChl 



molecules, the unit of B800 and B850 ring, are coupled 
uniformly to the incident light of frequency u) which is 
described by the creation (annihilation) operator (b) . 

To mimic the major function of natural process, 
we adopt the uniform couplings in the present model, 
namely, 6j = e. Then the model Hamiltonian reads as 
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H = H D +e A A f A+H DA +ujb f b+jy] (e+b + h.c) , (1) 
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\e) A (g\. The physical implementation of the 



present structure could be quantum dot, which has the 
size about 2 — lOnm and distance about lOnm [lj|. The 
recent experimental synthesis of 12-porphyrin ring [20| 
also opens up the possibilities of designing resemblance 
of natural light-harvesting element. The generic model 
with N donors and M acceptors has been discussed with 
the master equation Q . The similar setup with all donors 
and acceptors in a chain has also been discussed to reveal 
the optimal constitution of the two components |13[ . 

In this letter, we only consider the effective transfer- 
ring process and explore the advantage of the spatial 
configuration of the ring type. The collective excitation 
of the donor ring is described by the Fourier transfor- 
mation ej = Ek etkJ ek/VN and ej = £ fc e-**»'4/V^V, 

where eu {^Cj is the annihilation (creation) operator of 
the collective mode with definite momentum k. Here, 
the summation is over all the discrete momentum k n = 
2tt(ti — 1)/N, with n — 1,...,N. Indeed, in the large 



N— limit, we can show that [efc,e 



>kk' 



thus the 



collective excitations behave as bosons [21] . In terms of 
the boson- like operators eu \^h) > ^ ne donor Hamiltonian 

is seemly diagonalized as Hjj = ^ fc (e + 2g cos £;) e]],e~fc, 
which represents an energy band with N sub-energy- 
levels. We note that interaction term is rewritten as 
Ho a = V~Nto {cqA^ + h.c), which shows that only the 
zero mode excitation with k — is coupled to the ac- 
ceptor while the others described by fe]Q (k ^ 0) 

are decoupled with the acceptor. Therefore, the above 
Fourier transformation separates the total Hamiltonian 
into two un-coupled parts, H' = Hd — (e + 2g) ejeo and 
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where ujq = e + 2g — in and to a — — £T. We need to 
point out that the dissipation rates k and T have been 



phenomenologically introduced to describe the loss of ex- 
citations from the donors and the acceptor respectively. 

Next we give two remarks on the implication of the 
above effective Hamiltonian: 1. the zero mode excitation 
described by ej has the energy Ek=o = e + 2g. Acting on 
the ground state |0) = |<7i, #2, --jSiv); gi yes a uniform 
superposition |1) = |lfe=o) = Ej=i /V^V of the single 
localized excitations |lj) = \gi,...,gj—t,ej,gj+i,...,gN) 
in the jth donor (J — 1,2,..., N). It is similar to the sin- 
gle magnon state in the spin wave system. It has been 
numerically proved that the collective excited initial state 
|lfe=o) would result in the maximum efficiency with re- 
spect to any other mode [J]. 2. In single excitation sub- 
space, the photon-assisted donor-acceptor system could 
be described as a three-level A system as illustrated in 
FigQJb). The dash-boxed area is the corresponding sin- 
gle excitation subspace. Interestingly, the incident light 
only couples to the zero mode of the donor ring. Since 
all other modes are decoupled from the capture process, 
they do not contribute to the light-harvesting process, 
thus the excitation energy is transferred to the accep- 
tor only through the zero mode channel. Therefore, the 
transferring efficiency will be improved dramatically, if 
the nature light-harvesting systems were optimized to 
emerge a zero mode. 

In the single excitation case, the evolution of excita- 
tion is constrained in the subspace spanned by = 
|1,0,0),|1 D ) = |0, 1,0) and \l A ) = |0, 0, 1) for the direct 
product state \a, b, c) = \a)(£)\b)<g)\c) of the photon, donor 
in zero- mode and acceptor respectively. Let \(j>(t)) be the 
single excitation wave function with corresponding am- 
plitudes u(t), v (t) and w (t) to the above basis vectors. 
The Schrodinger equation is reduced into iV (t) = MV (t) 
for V (t) = [u (t) ,v(t),w {t)] T and 



M = 




(3) 



The energy transfer is usually understood as the de- 
cay from the donors in excited states 0, [Hj], thus the 
overall transfer efficiency is given by an integral rj = 
/ °° 2r \ v (t)\ 2 dt. In practice, with the hopping from the 
donor to the acceptor, the rate of transferring excitation 
to outside agent should be larger than that only by the 
dissipation from the excited donors, i.e., V^/Vio ^ k and 
r > k. Typically, we will choose the parameters 
here as e^/io = 10, g/to = 0.3 and T/to — 0.3 with 
to m lOps^ 1 . In most cases, we set the total donor num- 
ber N = 8. We illustrate the dependence of efficiency 
on detuning oj — e A and dissipation parameter k/T in 
Figj2ja). There exists a high peak at u = e A for given 
k/T with the peak value of efficiency almost unchanged. 
This observation reflects the efficiency is insensitive to 
the decay of the donor ring. At resonance, the exci- 
tation amplitude of the donor is highly suppressed, as 



3 



I' 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

K r 




FIG. 2: (Color online) (a) Transfer efficiency r\ vs detuning 
oj — €a and dissipation parameter n/T. The slow change of 
peak value indicates that the efficiency at resonance is insen- 
sitive to the donor decay rate k. (b) Amplitude evolution 
|u(t)|, \v (t)\ and \w (t)\, (c) Amplitude ratio \w(t)/u(t)\ at 
resonance uj = ea and k = T/3. 



illustrated in Fig|2]Jb). For the system comprises only 
the donor and the acceptor with initial excitation on the 
donor, there is a up-bound for the overall transfer effi- 
ciency jy^ax = r/ (r + k) . In the present discussion with 
the light capture included, the efficiency actually goes 
beyond the up-bound. 

The above discovery that the transfer efficiency is in- 
sensitive to the noisy of the donors at resonance can be 
explained according to the dark state, which has been 
widely investigated in quantum optics 22[. To this 
end, we write the dark state \D (t)) 



cos0(t)|l 6 ) - 

sin 6 (t) | 1a) in the single-excitation subspace for T = 0, 
where 8 (t) = arctan [J/to exp [— i (ca — u) t]]. \Dq (t)) 
is an eigenstate with vanishing eigen-value, namely, 
H e g \Dq (t)) = 0. It has been demonstrated that the 
perfect transfer of population between the two low-lying 
energy levels can be achieved by adiabatically tuning the 
Rabi frequencies J and to (23[. In this process, the ex- 
citation on the upper energy level is suppressed to avoid 
the dissipation of excitations. 

This dark stated based mechanism persists in the 
present artificial system with bionic structure. For this 
system, the evolution wave function is 



•(*)> = £ 
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\E, 



(4) 



1/2 [({Xi - LJ ) [Xi - u A ) - Nt 2 ) \l b ) - 

NJto | 1a)] is an eigenstate with a 



where \Ei) = Af 4 

\fNJ{Xi-UJ A )\lD 

normalized constant J\f 2 — N 2 J 2 t 2 + NJ 2 \x t - loa\ + 
\(xi — ujo) (xi — uja) ~ Nt^\ ; Xi is the corresponding 
eigenvalue of the matrix M, whose expressions are not 
explicitly written down since they are too lengthy. If we 
choose the parameters as previous e^/to = 10, g/to = 



0.3, K/t = 0.1, r/< = 0.3 and J/t = 0.1, the eigen 
states at resonance uj = ea are written as 

\E X ) ~ -0.995 |l 6 )+0.01i|li ? ) + 0.01 

\E 2 ) ~ 0.07 11 6 ) + (0.7 + 0.03i) 0-7 \1 A ) , 

\E 3 ) ~ 0.07|l 6 ) + (-0.7+0.03i)|lrj)+0.7|l A ).(5) 



Here, \Ei) is the dark state with very small component 
in donor excitation, which is proportional to JT/^/Ntl. 
The initial one photon state can be rewritten with the 
above eigenstates as the basis, i.e., |l) b = 0.995 \Ex) + 
0.07 \E 2 ) +0.07 \E 2 ). 

The component of dark state in the initial state is ap- 
proximate pE\ — W \Ao + J 2 : which is almost 1 un- 
der the practical condition J <C to- In the capture and 
transferring process, the population on excited donors is 
suppressed to be small in avoiding dissipation, as illus- 
trated in FiglJJb). On the dark state, the system de- 
cays very slowly at the rate j D ~ T (J/to) , while it 
dissipates quickly on the bright state (\E 2 ), \E 3 )) at rate 
7s ~ (k + T) /2. Thus, for large time scale tb > 1/^fB, 
the main contribution of transferring is carried on by 
the dark state. We demonstrate the ratio of populations 
on the photonic state and on acceptor in FigO^c). The 
asymptotic value of this ratio for long time is approxi- 
mately to I '\Ao + J 2 , which is the one of dark state in 
the initial state. 

In the transferring process, another important quan- 
tity is the average transfer time, which is defined as 
t = r; _1 J Q °° 2Tt \v (t)f dt. To effectively utilize the en- 
ergy, the excitation should be transferred with a high effi- 
ciency and also within a short time scale. In the previous 
discussions, we have proved that the efficiency can be im- 
proved via the dark state mechanism. However, we have 
illustrated that the decay rate of photon at resonance is 
suppressed by a factor (J/t ) 2 , while the efficiency is pro- 
portional to [1 + ( J/to) 2 ]^ 1 / 2 ■ We meet a dilemma that 
the efficiency and the average transfer time can not be 
optimized simultaneously. In FigJ3Ja), we demonstrate 
the average transfer time as a function of the dissipation 
rate of the donor and the detuning. The small peak in 
the center hints that the transfer time is not optimal at 
resonance uj — ca- It is readily seen in FiglH^a) that the 
optimal frequencies of quick transfer are not at resonance 
(lu = e + 2g) but split into two, which is known as the 
Rabi splitting in quantum optics. The two peaks can be 
determined by exactly diagonalizing the photon assisted 
donor-acceptor system in single excitation subspace as 
w± = (e A + 2g + e)/2±[(e + 2g-e A ) 2 /A + Nt 2 } 1 / 2 . The 
transfer time reaches its minimum optimal point when 
lu = u>±, while the efficiency is not at its maximum as 
illustrated in Figj3] 

In fact, such dilemma has been met in many investi- 
gations about heat engines the Carnot heat engine 
converts the heat into work with maximum efficiency, 
while it takes infinite long time. In practice, one would 
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FIG. 3: {Color online), (a) Average transfer time vs fre- 
quency detuning u> — ea- (b) mean power output vs de- 
tuning w — £a for different dissipation rate on the donor ring 
k = r, 1.5r and 2.or. 



concern more about the output power, which character- 
izes the output energy within unit time. In the present 
case, we introduce a similar quantity V = rj/r character- 
izing the ability of the excitation energy transfer, which is 
called the mean transfer power. In FigJ3](b), we demon- 
strate the mean power V as a function of dentuning be- 
tween the incident light frequency and the acceptor exci- 
tation energy for different donor dissipation rates k = T, 
1.5r and 2I\ The optimal frequency is uj = 9.1io and 
uj = 14.7io, which are different from the optimal value of 
both efficiency and average transfer time. Actually, we 
have observed that the photosynthesis systems choose to 
absorb sunlight mainly from two domains of the solar 
spectrum. For example, the green plants mainly use red 
and blue color photons. The present model can be uti- 
lized to justify this observation. However, for the practi- 
cal system, the estimation of the exactly frequency goes 
beyond the scope of the present model because it con- 
cerns very complicated biological environments. 

In summary, we have studied the light capture and 
excitation transfer process in a generic model consisted 
of donor and acceptor assisted by photons. By optimiz- 
ing this artificial photosynthesis system to realize an ef- 
fective three-level A configuration evolving dark state, 
we demonstrated the coherent population transferring 
through the dark state channels, where the dissipation 
from donors is effectively suppressed. In the present 
studies, we deal with the dissipation of the excitation by 
phenomenologically introducing an imaginary part to the 
Hamiltonian of the donor. However, in reality, this dis- 
sipation is always connected to the vibration degrees of 
freedom, which may be account for the dimerized struc- 
ture as discussed in the LHC I and LHC II [l5J. And also, 
it is worth to investigate effect of vibrations in designing 
some artificial light-harvesting systems. For the present 



model, we can also discuss the quantum or classical cor- 
relations of the output excitations with more than one 
photon, which could be used to explain the mechanism 
of the light-harvesting system avoiding the damage from 
high intensity light source. 
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